Echoes of the fifth dimension? 
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In this article we examine the question of whether the highest energy cosmic ray primaries could 
be ultra relativistic magnetic monopoles. The analysis is performed within the framework of large 
compact dimensions and TeV scale quantum gravity. Our study indicates that while this hypothesis 
must be regarded as highly speculative it cannot be ruled out with present data. 
PACS numbers: 98.70.Sa, 13.85.Tp, ll.lO.Kk 



O 
O 
O 

> 

O 

m 
> 

m 

o\ 
o 
o 
o 



Oh: 

CD , 



X 



This past year has seen a massive resurgence of inter- 
est in higher dimensional spacetimes [ll , a key new con- 
cept being the localization of matter, and even gravity 
on branes embedded in extra dimensions H- Depend- 
ing on the dimensionality and the particular form of this 
space, the long standing (Planck) hierarchy problem can 
find alternative solutions. In the canonical example of 
PI , the Planck scale of the four dimensional world is re- 
lated to that of a higher dimensional space-time simply 
by a volume factor. 
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where M, ~ 1 TeV is the fundamental scale of gravity, 
Mpi = 10^^ GeV, and n is the number of extra dimen- 
sions. With this factorizable geometry the case of one 
extra dimension is clearly excluded since gravity would 
then be modified at the scale of our solar system. How- 
ever, for n > 2, r is sufficiently small (the fundamental 
Planck scale is lowered all the way to the TeV scale) and 
the model is not excluded by short distance gravitational 
measurements. A more compelling scenario requires cur- 
vature to spill into the extra dimension B. Within this 
framework the background metric is not flat along the ex- 
tra coordinate, rather it is a slice of anti de Sitter space, 
due to a negative bulk cosmological constant balanced 
by the tension of two branes. In this non- factorizable 
geometry, the curved nature of the spacetime causes the 
physical scale on the two branes to be different, and expo- 
nentially suppressed in the negative tension brane. Such 
exponential suppression can then naturally explain why 
the physical scales observed are so much smaller than 
the Planck scale. Variants of this solution have been dis- 
cussed by many authors [pi . These models make dramatic 
predictions which can be directly confronted by current 
and future collider experiments [|6| , as well as cosmologi- 
cal observations W • The search for extra-dimension foot- 
prints in collider data has already started. However, as 
yet no observational evidence has been found |g]. 

Another seemingly different, but perhaps closely re- 
lated subject is the apparent lack of a high energy cut- 
off in the cosmic ray (CR) spectrum. Over the last few 
years, several giant air showers have been detected Q, 
with no sign of the expected Greisen-Zatsepin-Kuz'min 
(GZK) cutoff ^^. Initiated by single high energy par- 
ticles hitting the atmosphere, these are large pancake- 



shaped slabs of high energy particles which hit the ground 
at nearly the speed of light and can cover areas of many 
square kilometers. The origin and nature of the progen- 
itors is, at present, a deep mystery |12[ . Protons with 
energies above the GZK cutoff lose energy rapidly via in- 
elastic collisions with the cosmic microwave background 
(CMB) and thus presumably must come from a nearby 
source. This seems unlikely |l3]. A typical nucleus of 
the cosmic radiation is subject to photodisintegration 
from blue-shifted microwave photons, losing about 3-4 
nucleons per traveled Mpc ||l^ . Gamma rays of the ap- 
propriate energy have a short mean free path for creat- 
ing electron-positron pairs |14]. Although neutrinos can 
propagate through the CMB essentially uninhibited, at 
these energies the atmosphere is still transparent, and 
most of them interact in the Earth if at all. The dif- 
ficulties encountered in identifying a known particle as 
candidate have motivated suggestions in favor of "exotic" 
massive neutral hadrons, whose range is not limited by 
interactions with the CMB ||l^. However, the latter pre- 
dicts a correlation between primary arrival directions and 
the high redshift sources, which is not supported by the 
data set now available fl^. On a different track, it was 
recently put forward that extra dimensions may in prin- 
ciple hold the key to overcome this puzzle ||l7|. In this 
article we shall explore this fascinating possibility. 

It has long been known that any early universe phase 
transition occurring after infiation (say with symmetry 
breaking temperature T^), which leaves unbroken a t/(l) 
symmetry group, may produce magnetic monopoles p8[ . 
For instance, minimal SU{5) breaking may lead to "bary- 
onic monopoles" of mass M ^ Tc/a, with magnetic 
charge [/(1)em and chromomagnetic (or color-magnetic 
charge) SU {3)c 19 1. Here a stands for the fine struc- 
ture constant at scale Tc. These monopoles easily pick 
up energy from the magnetic fields permeating the uni- 
verse and can traverse unscathed through the primeval 
radiation. Thus, they are likely to generate extensive air 
showers pflj.n Before proceeding further, it is important 
to point out that if the monopoles are formed at the usual 



*The idea of monopoles as constitutents of primary cosmic 
radiation is actually quite old, it can be traced back at least 
as far as 1960 ^. 
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SURVIVAL OF THE GALACTIC MAGNETIC FIELD 



FIG. 1. The upper flux limit of monopoles (as a function 
of the relativistic parameter /3) as observed by the AMANDA 
experiment, is shown together with the Parker bound and 
the CR flux assuming that the primaries at the end of the 
spectrum {E > 5 x 10^^ eV) have masses ~ 10^" GeV. 



grand unification (GUT) scale ~ 10^^ GeV, the energy 
density overdoses the universe. Thus, to avoid this effect 
the symmetry breaking scale associated with the produc- 
tion of monopoles has to be shifted to lower energies. Re- 
markably, if the GUT scale is at ^ 10^ GeV, one would 
end up with an abundance of relativistic monopoles well 
below the closure limit, and yet potentially measurable 
to explain the tail of the CR-spectrum. In addition, for 
such a critical temperature the observed flux of ultra high 
energy CRs is below the flux allowed by the Parker limit 
p2[ . Moreover, the CR flux does not violate the upper 
bound for the monopole flux based on preliminary results 
quoted by the AMANDA Collaboration (see Fig. 1) [||. 
Unfortunately, contrary to the observed CR arrival direc- 
tions, the expected flux of relativistic monopoles is highly 
anisotropic, pointing towards the magnetic lines near the 
Earth Q. 

In the multidimensional models, the low-scale unifi- 
cation enables the production of light-mass monopoles, 
say M ^ 100 TeV. Furthermore, the physical embodi- 
ment of these theories allows a natural generalization of 
the 't Hooft-Polyakov monopole providing a convenient 
set of representations for Dl-branes ending on D3-branes, 
and consequently even lighter monopoles. Note, however, 
that direct searches at accelerators pretty much exclude 
masses below a few hundreds of GeV, whereas bounds 
stemming from quantum effects on current observables 
turn out to be ~ 1 TeV 123]. The light- mass monopoles 
could lose and gain energy as they random-walk towards 
the Earth. The maximum energy attainable before hit- 



ting the atmosphere is roughly 10^^ eV |26). Therefore, 
these "particles" would be ultra-relativistic, and the ex- 
pected flux has no imprint of correlation with the local 
magnetic field. 

To mimic a shower initiated by a proton the monopole 
must transfer nearly all of its energy to the atmo- 
spheric cascade in a very small distance. The large in- 
ertia of a massive monopole makes this impossible if 
the cross-section is typically strong, ^ 100 mb. Wick, 
Kephart, Weiler and Biermann (WKWB) [^ have re- 
cently pointed out that this problem can be avoided in 
models in which the baryonic monopole consists of q- 
monopoles confined by strings of chromomagnetic flux. 
To describe the interactions of such a monopole in air, 
WKWB have developed a model based on the four fol- 
lowing axioms: i) before hitting the atmosphere the 
monopole-nucleus cross section is roughly hadronic ao ^ 
■'^QCD (unstretched state), attaining a geometric growth 
after the impact; ii) in each interaction an 0(1) frac- 
tion of the exchanged energy goes into stretching the 
chromomagnetic strings of the monopole; iii) the chro- 
momagnetic strings (of tension T ~ Aq(-,j-,) can only be 
broken to create monopole-antimonopole pairs (a process 
highly supressed and consequently ignored); iv) the av- 
erage fraction of energy transferred to the shower in each 
interaction is soft AE/E = rj Ki Aqcd/A/- 

Generally speaking, in this set up the monopole will 
penetrate deeply into the atmosphere (the cross section 
is comparable to that of a high energy proton). How- 
ever, since the geometrical cross-section grows propor- 
tionally with the Lorentz factor 7, the interaction length 
(after the impact) shrinks to a small fraction of the depth 
of the first interaction. Stated mathematically, the un- 
stretched monopole's string length, L ^ A^^, increases 
by 5L = A.E/T. Recalling that nearly all of the ex- 
changed energy goes into stretching the color magnetic 
strings, the fractional increase in the length is bLjL = 7, 
yielding cti ^ (1 -I- 7)/Aq(-;j3. Now, the total mean free 
path after the iV-tli interaction reads. 
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where we have assumed a constant density of nucleons 
"■nuci ~ (4/3) TT ARq and we have used the approximation 
JN ~ (1 — AQCD/Af)^7 ~ 7- Here A stands for the 
mass number of an atmospheric nucleus, and Rq « 1.2 — 
1.5 fm. It should also be stressed that for N = -q^^ 
the approximation has an error bounded by limAr^oo(l — 
jy-i-jAf _ g-i^ Pqj. y^-i ^ -|^^ ^Yie total energy traveled 

between the first interaction and the r/^^-th interaction 
is then 
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Note that the mean free path for all secondary interac- 
tions is 0(1/7) compared to the first one. All in all, a 
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FIG. 2. Lateral distributions of charged particles and 
muons from AIRES simulations of a 100 EeV monopole with 
M = 100 TeV as a function of the distance to the shower core 
R. The error bars (obscured by the points themselves in most 
cases) indicate the RMS fluctuations of the means. 
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FIG. 3. Atmospheric cascade development of 300 EeV 
monopole induced showers, superimposed over the Fly's Eye 
data. The error bars in the simulated curves indicate the 
RMS fluctations of the means. 



baryonic monopole encountering the atmosphere will dif- 
fuse like a proton, producing a composite heavy-particle- 
like cascade after the first interaction. 

To examine the signature of such a cascade, we car- 
ried out a Monte Carlo simulation of monopole showers 
a la WKWB using the aires program (version 2.2.1) 
pTf . Specifically, several sets of proton "clumps" , each 
containing Af/AgcD, were injected at 100 km a.s.l with 
the first interaction point fixed according to the proton 
mean free path. The sample was distributed in the en- 
ergy range of 1 x 10^* eV up to 3 x 10^" eV, and was 
equally spread in the interval of 0° to 60° zenith angle 
at the top of the atmosphere. All shower particles with 
energies above the following thresholds were tracked: 750 
keV for gammas, 900 keV for electrons and positrons, 10 
MeV for muons, 60 MeV for mesons and 120 MeV for 
nucleons. The hadronic interaction was modelled with 
the SIBYLL package [p3 . The results of these simulations 
were processed with the help of the aires analysis pro- 
grams. 

The resulting lateral distributions from a vertically in- 
cident monopole of 100 EeV (7 = 10^) for muons and 
charged particles are presented in Fig. 2. A distinctive 
signature of this kind of shower is the great number of 
muons among all charged particles. This feature was ob- 
served in one not well understood "super-GZK" event 
p9[ . Roughly speaking, a magnetic monopole could then 
be a candidate primary for the highest energy Yakutsk 
event. However, WKWB-monopoles associated with a 
symmetry breaking at Tc ~ 1 TeV certainly cannot ex- 



plain all features of the data at the end of the spectrum. 
This is illustrated in Figs. 3 and 4. In Fig. 3 we show the 
longitudinal development of monopole showers superim- 
posed over the experimental data of the world's highest 
energy cosmic ray to date [ pO[ . To get some numerical es- 
timates we analyzed the data by means of a x^ test [Q . 
We assume that the set of measured values by Fly's Eye 
are uncorrelated (any depth measurement is independent 
of any other), and make use of the quantity 
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where q is the total number of points in the analysis, ax ■ 
is the error on the Xjth coordinate, Xj is the measured 
value of the coordinate, and aj the (hypothetical) true 
value of the coordinate. For masses of a few hundred 
GeV, the obtained x^ increases with rising mass from 
13.9 to 58.4. Our analysis indicates that masses above 
600 GeV are excluded at more than 99 % C.L. On the 
other hand, WKWB monopoles of masses around 200 
GeV become an alternative explanation for the Fly's Eye 
event. It is important to stress that a monopole mass 
« 200 GeV is not favored by D0 data |^^ , although one 
should keep in mind that these bounds are quite model 
dependent. Moreover, in view of the wide variety of un- 
certainties in the Fly's Eye event (the total error in the 
energy determination is 93 EeV |3Q]), one may be excused 
for reserving final judgment until more data is available. 
A better understanding of the present situation needs 
the analysis of the evolution of the shower maximum 
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FIG. 4. Average slant depth of maximum of showers initi- 
ated by monopoles. The error bars indicate the RMS fluctu- 
ations of the means. 



^max with energy. To this end, the charge multiplicity 
(essentially electrons and positrons) was used to deter- 
mine the number of particles and the location of Xmax by 
means of four parameter fits to the Gaisser-Hillas func- 
tion pTJ l . The situation is summarized by displaying the 
mean Xmax as a function of the logarithm of the pri- 
mary energy in Fig. 4. It is clear that despite its deep 
penetration, the monopole cascade develops much faster 
than a proton shower ||33|] . It can be seen by inspec- 
tion that the Xmax values produced by ultra-relativistic 
monopoles {E ^ 10^^ eV) with masses > 200 TeV are 
inconsistent with those reported by the Fly's Eye exper- 
iment, whereas the X^ax values of showers induced by 
lighter monopoles, M < 500 GeV, are within 1 standard 
deviation of the scarce "super-GZK" data |g[ . 

Whether or not the laws of physics should be for- 
mulated in more than four dimensions is still unclear. 
The possible existence of compact large extra dimensions 
brings with it low energy phase transitions p4| , providing 
a profitable arena for baryonic light-mass monopole pro- 
duction. If this is the case, the monopoles could be accel- 
erated to ultra-relativistic energies as they roam through 
space, inducing extensive air showers from time to time 
after hitting Earth atmosphere. As we have discussed in 
this article, the atmospheric cascade of monopoles with 
M ^ 500 GeV could reproduce quite well the main fea- 
tures of the recorded giant cosmic ray showers. Cer- 
tainly, more data is needed to test the WKWB hypoth- 
esis. Forthcoming ground arrays and satellites, such as 
the Auger Observatory [p5[, the next SCROD [^, and 
EUSO/OWL/AirWatch g7|-||], will help to increase the 
CR sample and more precise limits on the air shower ob- 



servables will be available, shedding light on the ideas 
discussed in this paper. 
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